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SIMMARY 


A method is presented for the analytical determination of the natural 
coupled frequencies and mode shapes of vibrations in the vertical plane of 
tandem helicopters. The coupled mode shapes and frequencies are then used 
to calculate the response of the helicopter to applied oscillating forces. 
Degrees of freedom included in the analysis are translation, pitching, and 
bending of the fuselage; translation, flapping, and bending of the blades; 
and translation of the engine. The method employs the Lagrange dynamical 
equations for free vibrations in conjunction with the kinetic and potential 
energies of the system in order to obtain the differential equations of 
motion for the coupled system which are then written in matrix form. The 
elements of the determinant of the matrix equation include the natural 
frequencies, mode shapes, and mass distributions of the uncoupled compo- 
nents of the helicopter (such as fuselage, blades, and engine) and permit 
the inclusion of experimental data of the uncoupled components in the 
evaluation of the coupled frequencies of the coupled system. 

The results of calculations made for a particular tandem helicopter 
show the variation of the coupled frequencies with rotor speed and indi- 
cate the changes in the coupled frequencies which result from changes 
in the structural properties of the individual helicopter components. 
Calculated response curves show how the vibrations of the cockpit vary 
with the frequency and method of application of oscillating forces applied 
at the hubs. Suggestions for further refinement of the analysis to aid 
in the determination of the coupled frequencies for more complex and real- 
istic systems are also presented. 


INTRODUCTION 


In many cases the development of the prototype of a particular heli- 
copter is handicapped because severe structural vibrations are encountered. 
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Such vibrations often limit the utilization of production articles because 
of fatigue of the component parts of the structure as well as of the opera- 
ting personnel. Available information indicates that the high levels of 
vibration are attributable primarily to two factors; namely, the high 
periodic content of the aerodynamic loading on the rotor blades and the 
amplification of the structural deformations due to proximity of resonance. 
Accordingly, it appears that the vibration levels may be reduced either by 
reducing the levels of the oscillating aerodynamic loads, by the inclusion 
of appropriate da mp ing, or by choosing the natural frequencies of the cou- 
pled structure so that they are as far as possible from the integral multi- 
ples of the rotor speed. 

The present paper relates to resonance amplification and, in partic- 
ular, to the determination of the natural frequencies, mode shapes, and 
response of the coupled helicopter structure as it exists in flight. Con- 
ventional methods such as those of references 1 and 2 may be used to pre- 
dict the natural frequencies of the various components of helicopters such 
as fuselages, blades, and engines. Reference 3, which reports on the ini- 
tial part of the present investigation, presents a method for the predic- 
tion of the natural frequencies and mode shapes of the structure which 
result when the various components are coupled together. The purpose of 
this paper is to elaborate on reference 3 and to present a method for 
utilizing the resulting coupled frequencies and mode shapes in computing 
the response of the helicopter to applied oscillating loads. 

The method is an energy method that involves chosen modes and fre- 
quencies of the uncoupled components which may be determined either by 
expressing the potential energy of the components in terms of given stiff- 
ness and mass distributions or, preferably, from experimental vibration 
tests of the components. The method is given in terms of a tandem- 
helicopter configuration but other configurations can be treated in an 
analogous manner. 

The present paper is concerned primarily with vibrations in the 
vertical plane - in particular, with fuselage translation, pitching, and 
bending; rotor blade translation, flapping, and bending; and engine trans- 
lation. The characteristic frequency determinant is derived, and calcu- 
lated values for the natural coupled frequencies and mode shapes for a 
particular tandem helicopter are presented. These mode shapes and fre- 
quencies are then used to dete rmine the vibration amplitudes of the cock- 
pit which occur when oscillating forces are applied at the front and rear 
rotor hubs. The results of some trend studies are also presented to show 
the effect of variations of the uncoupled frequency of the fuselage and 
of the Southwell coefficient of the rotor blades on the coupled frequencies 
of the helicopter. 
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SYMBOLS 


coefficients of front-rotor -blade modal deflections 
r-p 

— and x F , respectively, in. 

Rp. 


coefficients of frequency determinant (defined after 
eq. (11)) 


g 

coefficients of fuselage modal deflections 1, — , 

and y^, respectively, in. 

coefficients of rear-rotor-blade modal deflections 

— and Z]_, respectively, in. 

% 

structural stiffness of helicopter at station a, 
lb-in 2 

O 

structural stiffness of front rotor blades, lb- in. 
structural stiffness of rear rotor blades, lb-in. 2 

P 

structural stiffness of fuselage, lb-in. 

applied forcing function (see eq. (12)), lb 
amplitude of applied forcing function 

structural damping coefficient 

radius of gyration of mass of fuselage about front 
rotor hub, in. 

spring constant for engine mounts, lb/in. 

Southwell coefficients for blade first elastic flap- 
wise bending mode (defined after eq. (ll)) 

distance between front and rear rotors, in. 
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*cg 


*e 

L 

m 


nip 


m R 


Mp 

Mr 

Me 

Mf 

M f,l 


M r 


Pi q. 
Pr> Qr 


r F 


r R 


Rp 

r R 


distance from front rotor hub to center of gravity 
of fuselage, in. 

distance from front rotor hub to engine, in. 

overall length of helicopter, in. 

mass per unit length of helicopter at station a, 
lb-sec 2 /in . 2 

mass per unit length of front rotor blade, 
lb-sec 2 /in . 2 

mass per unit length of rear rotor blade, lb-sec 2 /in . 2 

mass per unit length of fuselage, lb-sec 2 /in . 2 

mass of front rotor blade, lb-sec 2 /in. 

mass of rear rotor blade, lb-sec 2 /in. 

mass of engine, lb-sec 2 /in. 

mass of fuselage, lb-sec 2 /in. 

effective mass of fuselage (defined after eq. (ll)), 
lb-sec 2 /in. 

effective mass of helicopter in rth mode (see 
eq. (17) ), lb-sec 2 /in. 

number of blades on front and rear rotor, respectively 

real and imaginary parts of response of helicopter 
in rth mode (see eqs. (23) and ( 24 )), in. 

radial position of any blade element on front rotor 
(see fig. l), in. 

radial position of any blade element on rear rotor 
(see fig. 1 ), in. 

radius of front rotor, in. 

radius of rear rotor, in. 
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longitudinal position of any fuselage element (see 
fig. 1), in. 

time, sec 

centrifugal force on blade at any radial station, 
lb; kinetic energy, lb-in. 

centrifugal force on front and rear rotor blades, 
respectively, lb 

centrifugal force divided by square of rotor speed, 
lb-sec2 

potential energy, lb -in. 

deflection of engine, in. 

deflection of engine with respect to fuselage, in. 

deflection of element of front rotor blade, in. 

.uncoupled first-bending-mode shape for front rotor 
blade 

real and imaginary parts of response of helicopter 
(see eq. ( 26 )), in. 

deflection of element of fuselage, in. 

uncoupled first-bending-mode shape for fuselage 

deflection of element of rear rotor blade, in. 

uncoupled first-bending-mode shape for rear rotor 
blade 

coefficient of deflection of rth coupled mode of 
helicopter (see eqs. ( 13 ) and ( 23 )), in. 

amplitude of o^, in. 

distance measured rearward from leading edge of front 
rotor disk (see fig. 1), in. 

value of a where oscillating force is applied 



6 

NACA TN 3849 

0 

deflection of helicopter in response to applied 
force F (see eqs. (13) and (25)), in. 


rth-coupled-mode shape of helicopter 

M 

amplitude of 0 (see eq. ( 26 )), in. 

(0 

natural frequency of coupled modes; frequency of 
applied oscillating force, radians/sec 

“F,l* %1 

first natural uncoupled flapwise bending frequency 
of front and rear rotor blades, respectively, 
radians/sec 

“f,l 

first natural uncoupled vertical bending frequency 
of fuselage, radians/sec 


natural vertical frequency of engine on its mounts, 
radians/sec 

0* 

natural frequency of rth coupled mode of helicopter, 
radians /sec 

Q, 

rotor speed, radians /sec 


The deflection in the first bending mode of the fuselage, measured* 
at station s , is denoted by y^(s) . For example, y-]_(0) is the deflec- 
tion of the fuselage in first mode bending at the front rotor hub. (See 
fig. 1.) 


ANALYSIS 


Equations for Free Vibrations 

General considerations .- The analysis presented in this portion of 
the paper deals with the determination of the undamped, natural coupled 
frequencies of the vertical vibrations of a tandem helicopter. The fol- 
lowing degrees of freedom are considered: vertical translation, longi- 

tudinal pitching, and structural bending of the fuselage; vertical transla- 
tion, flapping, and flapwise bending of the rotor blades on the front and 
rear rotors; and vertical translation of the engine. This analysis does 
not treat fuselage side bending and torsion; however, an analysis of the 
structure for side bending and torsion can be made in a manner parallel 
to the present analysis. Although the analysis is given in terms of a 
tandem helicopter configuration, other configurations may be treated in 
an analogous manner . 
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The method follows the Lagrange energy-equation approach wherein 
the derivation of the equations of motion consists of writing the kinetic 
and potential energies of the structure in terms of the displacements w, 
x, y, and z. (See fig. 1.) These displacements are in terms of gen- 
eralized coordinates and chosen uncoupled modes of the structural com- 
ponents (fuselage, rotor blades, and engine). The energy equations are 
substituted in Lagrange's dynamical, equations for free vibrations to 
obtain the equations of motion which yield the coupled frequencies and 
mode shapes. 

Energy equations .- The kinetic energy T and the potential energy V 
are, respectively. 



( 2 ) 


Special attention is called to the fact that the energies of the front 
rotor with p-blades may be different from those of the rear rotor with 
q-blades. For this reason, in the derivation of the energies, the deflec- 
tions of the front and rear rotor blades are denoted independently by x 
and z, respectively. It is assumed in the present analysis that all the 
blades on a given rotor behave in a similar manner; however, equations (l) 
and (2) are readily extended to treat each blade as an independent degree 
of freedom. 

Choice of modes .- The next step toward obtaining the differential 
equations of motion is to express the deflections w, x, y, and z in 
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terms of chosen mode shapes. The choice of the translation mode, the 
pitching mode, and one bending mode for the fuselage; the flapping mode 
and one flap-wise bending mode for the blades on each rotor; and one 
translation mode for the engine leads to the following equations for the 
deflections : 

For engine translation: 


w = b 0 + b 0 y + b iyi (l e ) + w f 


(3) 


For front -rotor -blade translation, flapping, and bending: 

x = b 0 + b iyi (0) + — + a lXl (4) 


For fuselage translation, pitching, and bending: 


y = bQ + bQ j + bqyi 
L 


(5) 


For rear -rotor-blade translation, flapping, and bending: 

z = b 0 + b 0 + b iyi (l) + C 0 -- + (6) 

r R 


where the coefficients sq, aq, bQ, bQ, b]_, cq, Cq, and w^ are 

unknown functions of time, and x-[_, y -j_, and Zj_ are the chosen spatial 

functions or mode shapes. If higher modes are desired, it is only nec- 
essary to add appropriate terms to equations (3) to (6). 

Substitution of equations (3) to (6) into equations (l) and (2) leads 
to the following equations for the kinetic and potential energies, 
respectively: 
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-12 


b o + b iyi(°) + a© — + % X 1 
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ax F + 
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b o + b o + b iyi(^) + c o ~ + °i z i 


Rr 


dr. 


R 


(7) 


and 


V = i f * (EX) f (b iyi ") 2 ds + i Ke(v f ) 2 + J f ^ 


'0 


a l x l )^ dr ]p + 


^R O pRji / g- \ 

( e I)r(c 1 z 1 ") dr R + | + e-lX-l’] dr p + 


^R / c 0 \ 

i J 0 MS**!*!')** 


q 


(8) 


■where the dots denote derivatives with respect to time and the primes 
denote derivatives with respect to space or length. 

After the energies have been calculated, the next step is to apply 
Lagrange's equation for free vibration to obtain the differential equa- 
tions of motion; that is, 

±(*L) - + £v_ = o 

dt\8ao ) daQ 

l 

d /8 t \ _ + 8V_ = 0 

dtySa^i 8a^ 8a^ 


(9) 
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In the application of Lagrange’s equation (eqs. {9)), great simpli- 
fication is obtained if the chosen modes axe the natural modes of the 
uncoupled components; the rest of this paper is based on this choice. 
Equations (7) and (8) are then substituted into equations (9) to obtain 
the differential equations for free vibration. The characteristic equa- 
tions for free vibrations are found by considering harmonic motion, that 
is. 


aQ = “Sq sin cut 
a l = a l S ^ n f 


( 10 ) 


j 


where the amplitudes of a q, a^, . . . are denoted by a q, aj_, . . . . 

The result of substituting equations (7) and ( 8 ) into equations (9) and 
maVing use of equations (10) leads to the characteristic equations which 
are given in matrix form as follows: 


M f + Mg + 

pMp + qMR 

Mg -^ + qMg + 

% if 

AO 

Bo 

MeY^ie) 

iMjpy^O) 

^syiU) 

+ 

+ 

Al 


®1 

'l e 

Mg-— + qM R + 


0 

% 

Me^y^le) + 

0 


B 1 


+ <3Mr 








*0 

q 

4-(») 2 ] 

0 

y x (o)Ao 


0 


0 

Bq 

Eo 

0 


j yiU)Bo 


0 


0 

pMjy^O) + 

qMRyi(l) * 

Meyi( ? e) 

Me ~ + 

We? 1 ( 1 ) 

y^oJAo 

yi( l )Bo 

M e y 1 2 de) + 
pMpy 1 2 (o) + 
qM R yi 2 U) 

1 

1%,! + 

yi<0)Ai 


yjUJBi 

Al 

0 

0 

0 

yx(o)Ai 


4 -( ¥) 

_ 

0 







KF ' 1 ©i 



B 1 

Bl 

0 

0 

y^Bi 


0 

D 1 

t-M 









4)1 

Mg 

l e 

Me T 

0 

0 

MeVxde) 


0 


C 


Me 


b Q 



*0 

0 


*0 

0 



M^Ue) 


bl 

! 

0 


*1 

0 


*1 

"■ j 1 - (w)_ 


-f 


(n) 
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where 


Ao 




dr-r 


q r R R 

"o'Ei 0 “ RrR ^ 


*1 



mpXi dr p 


»l 



m R z l 


dr R 



D 0 




m R r R 2(Jr R 


C 


1 



m F x l 


2 dr p 


% 



m R z l 2dr R 


K F,1 


P_ 

Cl 



Tf( x 1* ) ^F 


K R,1 “ 


JL 

% 



M. 





nipy^ds 



Ke 

Me 


In the abbreviated notation used in the appendix, equation (ll) can be 
written as 



(11a) 


Coupled Frequencies and Mode Shapes 

The values of cu, the natural frequencies of the coupled modes of 
the helicopter, are determined by setting the determinant of equation (ll) 
equal to zero. The coupled mode shapes are determined by substitution 

of the modal coefficients t>o, Sq, and so forth, obtained from equa- 
tion (ll), into equations (3) to (6). 
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Inspection of equation (ll) reveals that, although there are eight 
degrees of freedom, the frequency term o> appears in only six of the 
diagonal terms when ft ^ 0 and in only four of the diagonal terms when 
ft = 0. Thus, the amount of work involved in the solution of the deter- 
minant can be reduced by changing the order of the matrix from eight to 
six when ft ^ 0 and from eight to four when ft = 0. The steps involved 
in this transformation as well as the final matrices obtained for ft = 30 
radians per second and ft = 0 for the basic configuration are given in 
the appendix. The appendix also contains some pertinent remarks regarding 
methods for obtaining the modal coefficients. 


Equations for Response of Helicopter to Applied Loads 

Once the natural frequencies and shapes of the coupled modes are 
obtained, the next objective is to determine the response of the heli- 
copter to applied aerodynamic loads. The calculation of the response 
necessitates a consideration of the damping of the structure . In the 
analysis that follows, a method is presented which permits the coupled 
mode shapes and frequencies determined by the methods presented in pre- 
vious sections of this paper to be utilized, with the inclusion of struc- 
tural damping, in the determination of the response of the helicopter to 
arbitrary forces. 

If $(d,t) be designated as the deflection of the helicopter at 
station a at time t (fig. l) , the differential equation of motion 
for an element of the helicopter at station <7 is 


(l + igh)- — yEI — — + m0 — F(d,t) (12) 

bo 2 \ so 2 / 


where F is the applied force, g h is the conventional structural 

damping coefficient, and simple harmonic motion is implied. If it is 
assumed that the total deflection can be expressed as a superposition of 
the natural coupled modes of the structure, then 


r=n 

0(cr,t) = 21 a r(‘t)0 r (cr) 

r=l 


(13) 


where 0 r (cr) is the rth- coupled-mode shape and a^t) is the time- 

dependent coefficient of the rth mode. Substitution of equation (13) 
into equation (12) yields 
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-®L(t)ia^i - a2(t)m^2 - ... - c^(t)n$ n + F(a,t) 




Since the natural modes of vibration are defined by 


ai-L tk) 

bA be 2 I 


= 


( 15 ) 


equation (l4) can be written as 


+ 1 Sh^ a l( t ) a b 2m ^l + + ig^ a 2(' t ) a h 2m ^2 + • • • + + 


- a2(t)m^ 2 - ... - ok(t)n$ n + F(ff,t) 


(16) 


Upon multiplication of equation (l6) by <f v , integrating over the length, 
and observing the conditions of orthogonality for natural modes, namely. 



equation (l6) becomes 


( 17 ) 




O « • 

ctro^. M r = -Oj-Mj. + 



F(cr,t)0 r dcr 


( 18 ) 


In general, F will be a distributed force and, if necessary, the 
integral term of equation (l8) can be evaluated. For purposes of illustra- 
tion, however, it is assumed that F is a concentrated force and the rest 
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of the analysis is based on that assumption. The integral term of equa- 
tion (18) is 



d a = F0 r (cr A ) 


( 19 ) 


where F is the force applied and 0 r ( cr A ) is the deflection of the rth 
mode at the point of application of the force. Therefore, 


^1 + igh^Or^^r + %% = F 0r( a A) 


( 20 ) 


If ccr and F are expressed as the real parts of a r ,o e ^ aJ k Fqs , 

respectively, then 


^1 + ig] a ^ a r,0 a r^Mr - a^ar^oMr - FQ^r(^A) 


The amplitude of is 


( 21 ) 


^r(^A) 


p 

(%• My 



( 22 ) 


and, by definition, the coefficient of the motion in the rth mode is 
given by the real part of the following expression: 


F (A-( a A) 


ac^Mj. 



+ iSh 


e iajfc = 


(p r - iQ r ) e ia * 


CXry — 


( 23 ) 
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where 


P r = 



Qr = 


Fo^r^A) 


o* 2 M r 



A 



(24) 


Substitution of equation (23) into equation (13) yields the total 
response 0 

the amplitude of which is 






+ 





(26) 


Then | 0 | is the amplitude of vibration at any station a in response 
to a harmonic force, applied at a = o^, which has an amplitude Fq and 

a frequency u>. ■ It should be noted in the application of equation (26) • 
that $ 2 _, $ 2 .’ • • • 0n are the values of the coupled modes at station o. 

In the event that it is desired that the response at any station be 
calculated for loads applied at several different stations, then some 
reduction of the work required may be obtained by suitably modifying the 
coefficients P r and Q r . Assume, for example, that (p r ) and (Qr)j 

are the in-phase and out-of -phase components of the response in the rth 
mode to an oscillating load applied at the front hub and designated by 
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subscript I. From equation ( 26) , the response in the rth mode at the 
same station to a load of the same magnitude applied in the same direc- 
tion at the rear hub, and designated by subscript II, is given by the 
components : 


( p r ) jj - ( p r ) -fix (<*a) jj jtfA o A ) j 


(27) 


The in-phase and out-of -phase components of the response due to all the 
modes is then given by 


r=n 

X II = Z(P r ) *r 
r=l II 

> 

r=n 

Y II := 2Z (Qr)n^r 
r=l 


(28) 


The response due to any combination of forces applied at either or both 
hubs in the same or opposite directions is then obtained by adding Xj, 

Xji, Yj, and Yjj in the proper sense and proportions. 


DISCUSSION OF RESULTS 


Calculation of Coupled Frequencies and Mode Shapes 

Scope of the calculations . - The analysis derived in this paper was 
used to calculate the natural coupled frequencies of a tandem helicopter 
having the structural parameters given in table I. The resulting modal 
coefficients are given in table II and the natural frequencies and mode 
shapes are given in figures 2 and 3* Some additional calculations were 
made to evaluate the effects of the natural frequency of the uncoupled 
fuselage system and the Southwell coefficient of the blades on the natural 
frequencies of the coupled system. These results are shown in figures 4 
and 5. Calculations were also made for a three -degree -of ^ -freedom system 
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(bending of the fuselage and bending of the blades on the front and rear 
rotors) and the resulting frequencies are compared with appropriate fre- 
quencies for the eight-degree-of -freedom system in figure 6 to obtain 
some idea of the relative contributions of the various degrees of freedom 
to the overall problem. 

Frequencies and mode shapes for the basic configuration .- The basic 
configuration is designated by the parameters given in table I which 
correspond closely to those of an existing tandem helicopter. The natural 
frequencies and mode shapes of the coupled system were calculated for 
this configuration for three rotor speeds: namely, 0, 20, and 30 radians 

per second. The natural frequencies are plotted in figure 2 as a func- 
tion of rotor speed, and the modal coefficients are tabulated in table II 
for each of the six natural frequencies and for each of the three rotor 
speeds. A rotor speed of 30 radians per second corresponds very closely 
to the normal operating -rotor speed for the chosen helicopter. The 
modal coefficients presented in table II for this speed were used in equa- 
tions (3) to (6) to obtain the mode shapes which are presented in figure 3 
and which are also sketched in figures 2, 4, and 5 to identify the modes 
of vibration that correspond to the different natural frequencies of the 
coupled system. 

The effect of coupling of the various degrees of freedom on the 
natural frequencies is indicated in figure 2 by a comparison of the 
coupled frequencies (shown by solid lines) with the uncoupled frequencies 
(shown by short-dashed lines). The curve for uncoupled blade first bending 
represents the blades of both the front and rear rotors. There are also 
sho.rt-dashed curves for uncoupled blade flapping; however, these curves 
cannot be seen inasmuch as they are coincident with the frequency curve 
for the lower coupled mode. The frequencies of the two- and three-per- 
revolution harmonic component of the aerodynamic loading (designated 2-P 
and 3-P herein and indicated by the long-dashed lines in fig. 2) are 
also presented to permit comparison with the natural frequencies. With 
the exception of the engine translational frequency and the fuselage 
first-bending frequency, the differences between the coupled and uncoupled 
frequencies are negligible at low rotor speeds. As the rotor speed is 
increased, the uncoupled frequencies of the fuselage mode and the blade 
bending modes approach each other, and, as expected, the effects of 
coupling become pronounced. 

The helicopter used as an example for these calculations has three 
blades on each rotor, and, therefore, the 3-P components of the aero- 
dynamic loading on the rotor blades are additive at the respective hubs. 
Thus, the region of primary importance on the frequency diagram of fig- 
ure 2 is the region at or near the crosshatched marks on the 3-P line 
which brackets the normal variation of the design rotor speed. 
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Inasmuch as the two rotors rotate in opposite directions with a 
phase difference of blade position on the two rotors of approximately 6o°, 
blades from both rotors are advancing into the wind simultaneously three 
times per revolution of each rotor. This condition suggests that the 
resultant 3-P periodic forces at the two hubs are in phase and would 
be particularly effective in the excitation of the symmetrical modes of 
the helicopter which have natural frequencies in the general region of 
85 radians per second. Figure 2 shows that there are three coupled modes 
of the helicopter which have natural frequencies reasonably close to the 
3-P frequency, two of which are symmetrical and the other antisymmetrical. 
In view of the preceding discussion, the data of figure 2 indicate that 
at the higher rotor speed (30 radians per second) the response of this 
helicopter in the third symmetrical mode to the 3-P oscillating excita- 
tion forces would be increased substantially by resonance amplification. 
Furthermore, if there is a component of the 3-P aerodynamic loading at 
the two rotor hubs which is out of phase, it appears likely that the 
structural response in the second antisymmetrical mode would also be 
increased by resonance amplification. 

The mode shapes presented in figure 3 show that the mode shapes 
alternate between symmetrical and antisymmetrical configurations . The 
relative motion of the engine with respect to the fuselage is shown to 
be negligible for the lower frequency modes but becomes relatively large 
at the highest frequency mode (sixth- coupled-mode shape). 


Effect of Variations in Uncoupled Components 

Effect of natural uncoupled frequency of the fuselage .- The analysis 
presented in this paper permits the inclusion of the natural uncoupled 
frequency of the fuselage as one of the parameters. This frequency can 
be calculated by analytical methods such as those of references 1 and 2; 
however, this requires an accurate knowledge of the mass and stiffness 
distribution of the fuselage which may be difficult to obtain. The accu- 
racy to which this frequency can be obtained experimentally for a given 
helicopter depends upon such considerations as the manner of support and 
the type of shaker installation. Since this frequency is one of the 
primary parameters, it is of interest to determine to what extent the 
natural frequencies of the coupled modes are changed by variations of 
the uncoupled frequency of the fuselage. Conversely, if the natural 
frequencies of the coupled system for a known design value of the natural 
uncoupled frequency of the fuselage are of such a magnitude as to lead to 
resonance amplification of the structural response, it is desirable to 
evaluate the extent to which the coupled frequencies can be changed by 
varying the uncoupled frequency of the fuselage. 

An indication of the change in the natural frequencies of the coupled 
modes due to changes in the uncoupled frequency of the fuselage is shown 
by the curves of figure 4. The calculations show that the frequencies of 
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the second and third symmetrical modes are changed appreciably whereas 
the frequencies of the antisymmetrical modes are not affected- Since 
the amount of fuselage bending in the antisymmetrical modes is negli- 
gible in comparison with the amount of fuselage bending in the symmetrical 
modes, the results shown in figure 4 are in agreement with expectations. 

Effect of Southwell coefficient for the blade first elastic flapwise 
bending mode .- The value of the Southwell coefficient chosen for the blade 
first elastic flapwise bending mode has a direct effect on the uncoupled 
frequency of the rotating blade. An indication of the effect of this 
parameter on the coupled frequencies is shown in figure 5* The value 
of Kp q estimated by means of reference 4 for the blade considered 

is 6.2 and the range of values chosen for these calculations is believed 
to be more than adequate to allow for estimation errors. These calcula- 
tions show that the coupled frequencies of both the symmetrical and anti- 
symmetrical modes which have frequencies approximately equal to the uncou- 
pled blade-bending frequencies are increased slightly as the Southwell 
coefficient is increased, and that errors in the values of the coupled 
frequencies due to the estimation of the Southwell coefficient Kp^q 

should not exceed 2 or 3 percent. 


Comparison of Results of Three-Mode and Eight-Mode Analysis 

Figure 2 shows that the uncoupled modes which are primarily affected 
by coupling in the frequency region of the 3-P excitation forces are 
fuselage bending and the bending of the blades on the front and rear 
rotors. Calculations were made to determine how well the frequencies of 
the coupled modes could be predicted by using only these three uncoupled 
modes. The calculation procedure consisted of retaining only the fifth, 
sixth, and seventh rows and columns of the frequency determinant of equa- 
tion (ll). The results are given in figure 6 together with the corre- 
sponding frequencies obtained from the eight-mode analysis. A comparison 
of the frequencies for the two cases shows that the frequencies calculated 
by the three-mode analysis are from 3 to 6 percent lower than those 
obtained by the eight-mode analysis, and suggests that a three-mode anal- 
ysis might be useful for obtaining a first approximation of the coupled 
frequencies in this region for configurations similar to the one treated 
herein. 


Calculation of the Response of the Helicopter to Applied Loads 

The methods developed in the previous sections of this paper were 
applied to calculate the amplitudes of the response of the helicopter 
at the location of the cockpit for harmonic oscillating forces applied 
at the front and rear rotor hubs. The amplitudes of vibration were cal- 
culated as a function of frequency of the periodic force for one value 
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of damping (gh = 0.1) and at a rotor speed of 30 radians per second for 
the four cases which follow: 

Case I: Oscillating force of unit magnitude applied in a vertical 

direction at the front rotor hub. 

Case II: Oscillating force of unit magnitude applied in a vertical 

direction at the rear rotor hub. 

Case III: Oscillating force of one-half unit magnitude applied in 

a vertical direction at both the front and rear rotor hubs in the same 
direction. 

Case IV: Oscillating force of one-half unit magnitude applied in 

a vertical direction at the front and rear rotor hubs in opposite direc- 
tions (producing a pure pitching moment on the fuselage) . 

The in-phase and out-of -phase components of the response, or cock- 
pit deflections, are given in figures 7 and 8 as a function of frequency 
for cases I and II, respectively. The natural coupled frequencies are 
indicated by the vertical dashed lines at the bottom of the figures. 

The response curves are presented to familiarize the reader with the 
nature of the components and as an aid in the calculation of the cockpit 
deflections due to other combinations of loads at the front and rear rotor 
hubs. The cockpit deflections for cases III and IV as well as for cases I 
and II were obtained by suitably combining the in-phase and out-of -phase 
components of figures 7 and 8 as outlined in the analysis section. 

The amplitude of the cockpit vibrations resulting from the oscillating 
forces applied to the hub in the four cases previously outlined are given 
in figures 9 to 12, respectively. Before discussing these figures, the 
reader's attention is again directed to figure 2 where, at the rotor speed 
of 30 radians per second for which these response calculations are appro- 
priate, six distinct natural frequencies are indicated within the range 
from 0 to 180 radians per second. These natural frequencies are indi- 
cated by the vertical dashed lines at the bottom of figures 9 to 12. For 
the hypothetical case of zero damping, it is expected then that the 
response curves similar to those of figures 9 to 12 would indicate six 
peaks of infinite amplitude at frequencies corresponding to the natural 
frequencies. With the inclusion of a small amount of damping, these 
peaks become finite and as the damping is progressively increased some 
of the peaks disappear. The peaks which remain when the forces are 
applied in a given manner indicate the modes in which the given forces 
are more effective in doing work. 

The curves of figures 9 to 12 show the responses obtained for = 0.1 
for the four cases investigated. The manner of application of the load in 
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each case is indicated in the figure by the accompanying sketch of a 
tandem helicopter. In a few instances where the number of calculated 
points, indicated by the small circles, are insufficient to define the 
peaks completely, the peaks are indicated by the long-dashed lines. All 
the figures indicate peaks at frequencies of approximately 32, 72, and 
92 radians per second and show that the coupled modes which contribute 
most to the responses are the second, third, and fifth coupled modes, 
respectively, of figure 3. This effect was also observed during the cal- 
culations by inspection of the values of the individual values of P r 

and Q r . (See eq. (24)). The most significant deviation from the gen- 
eral trend of high responses at forcing frequencies near the natural fre- 
quencies of the aforementioned coupled modes is found in figure 12. In 
this case, the response to forces applied so as to produce a moment on 
the helicopter is found to diminish considerably at the higher forcing 
frequencies. This reduction of the response at the higher frequencies 
appears significant because it indicates that, if the aerodynamic forcing 
functions which are usually encountered in flight could be forced to 
occur in an antisymmetric manner, the response at a frequency of 3-P 
or 90 radians per second would be substantially reduced over that obtained 
for symmetrically applied loadings as shown by figure 11. 


Perhaps one of the most significant points brought out by these 
studies is the fact that it is not always possible to obtain the natural 
frequencies of all the coupled modes of a structure by observing the 
response of the structure to arbitrary inputs. Conversely, if the applied 
loads are realistically chosen to be representative of those encountered 
under normal flight conditions, the peaks of the response curve appear to 
give fairly good indications of the natural frequencies of the more impor- 
tant coupled modes. 


Considerations Regarding Further Refinement of the Method 

The studies presented in this paper - although believed to be useful 
in the evaluation and reduction of the critical vibrations of tandem heli- 
copters by emphasizing the importance of coupled natural frequencies and 
mode shapes, providing a method for estimating them, and, in turn, using 
them to determine the response of the helicopter to applied loads - do 
not provide a complete answer to the problem. For one thing, the number 
of degrees of freedom treated, although they are believed to be the more 
important ones, are not nearly sufficient to predict all the coupled fre- 
quencies. A more complete analysis would also include the fore-and-aft 
motions of the front and rear rotor hubs, the pitching motions of the 
engine, and the chordwise motions of the blades in bending and rotation 
about the drag hinge. In fact, each blade should be treated as an inde- 
pendent flexible structure, subject only to the boundary conditions 
imposed by the vibrating hub to which it is attached. The complexities 
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of treating such a system analytically are apparent and they exceed the 
intent of the present paper. 

The method of analysis presented in this paper for the longitudinal 
vibrations - namely, the coupling of the structural properties of the 
uncoupled components of the helicopter to obtain the natural frequencies, 
mode shapes, and response of the coupled system - also appears applicable 
to the treatment of the case of coupled side bending and torsion of the 
fuselage. Vibrations of this type have been observed in tandem helicopters 
and a study of the effects of the significant parameters involved would 
provide useful design information. 


CONCLUDING REMARKS 


A method is presented for the calculation of the coupled frequencies 
and mode shapes of vibrations in the longitudinal plane of symmetry of 
tandem helicopters. The coupled frequencies and mode shapes are then 
used to calculate the response of the helicopter to. applied loads. Sample 
calculations were made for an existing helicopter and these calculations 
yield the coupled frequencies and mode shapes and the deflection of the 
cockpit to oscillating loads applied at the front and rear rotor hubs. 

The mode shapes, both symmetrical and antisymmetrical, show the relative 
deflections of the various components of the structure when excited in 
a given mode. 

• The results of the calculations for the helicopter chosen show that 
the natural frequencies for the coupled modes may differ considerably 
from the natural frequencies for the uncoupled modes; therefore, there 
may be a need for making a rather comprehensive coupled-mode analysis 
so that regions of adverse dynamic response can be detected and avoided. 

The results of calculations for different values of the uncoupled 
natural frequency of the fuselage indicate that the natural frequencies 
of the symmetrical coupled modes in the i general range of the uncoupled 
frequency of the fuselage vary appreciably with changes in the uncoupled 
frequency of the fuselage; however, the effect on the natural frequency 
of the antisymmetrical coupled modes is negligible. An increase in the 
values of the Southwell coefficient for the blade first elastic flapwise 
bending mode resulted in a small increase in the coupled frequencies of 
both the antisymmetrical and symmetrical coupled modes having frequencies 
in the proximity of the blade frequencies. 

The response calculations show that the amplitude of the vertical 
vibrations of the cockpit changes appreciably with changes in frequency, 
point of application, and manner of application of applied oscillating 
forces. High amplitudes occur when the frequencies of the applied oscil- 
lating forces are approximately equal to the natural frequencies of the 
second, third, and fifth coupled modes of the helicopter and the results 
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of calculations indicate that the contributions of other modes to the 
vibrations of the cockpit are secondary for the specific forces chosen. 

The results indicate that the vibration of the cockpit due to oscillating 
forces having a frequency of 3 per revolution (3-P) or 90 radians per second 
are substantially greater for forces applied symmetrically at the front 
and rear rotor hubs than for forces applied antisymmetrically. The results 
also emphasize that it is very unlikely that all natural frequencies of a 
helicopter could be obtained by observing the response of the helicopter 
to oscillating forces applied at any given point. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , August 2, 1956. 



24 


NACA TN 384-9 


APPENDIX 


STEPS IN REDUCING THE ORDER OF THE MATRIX EQUATION 


Procedure 'When Q ^ 0 

When Si 4 0, the matrix of the differential equations of motion, 
equation (ll) of the text, can be written in concise notation of sub- 
matrices as follows: 



'[»] M 

W M 


{’ll) 



(Al) 


where 


Mf + M e + pMp + qM R 


M e — + qM R + Mf ^ 

It It 


Me j+ qM R + Mf 
Mf(|f + Me/y) + qM R 


Ao Bo M e y 1 (l e ) + pMj-y-^O) + qM^Cl) Aq Bq Me 

[pj = (A3) 

0 Bo Mg j- yq(l e ) + lM R y -]_(*) 0 Bq Me y 
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j pj 1 = Transpose of J 


yi(o)Ao 


yi( l )Bo 


y 1 (0)A 0 y x ( 1 )Bq Uf fl + Mey^Cle) + y 1 (o)A 1 y 1 (l)B 1 Mey^lg) 

pMpy^lO) + 

< i%y 1 2 ( i ) 


y x (o)A 1 


yiU)®! 


^eyi(^e) 




Ci^.i) 2 + %,ia 2 ] 


D lj^(“R,l) 2 + K R,1 Q2 ] 



2 6 
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^1 = 


b 0 

/v 

bo 


■ 47 ) 


and 


m 


= < 


*0 

c 0 

^Vi 

b l 

a l 

ro 

C 1 

W f 


(A8) 


If [a] s a, j^} = rii, and so forth, then the algebraic expressions 

for the matrix equation, equation (Al) , are as follows: 


ai^ + 3^2 “ 0 


(A9) 


p'q-j. + bT l2 = ~~Z €T 12 
or 


(A10) 


If equation (A9) is solved for and the result substituted into equa- 

tion (A10), then 


( 6 




(All) 


Multiplication of equation (All) by the inverse of e reduces it 
to the desired form; however, the resulting matrix is unsymmetrical and 
poorly conditioned. These characteristics are undesirable if iteration 
procedures are used to determine the frequencies and mode shapes, parti- 
cularly so in this case because the frequencies are in some instances 
nearly equal to each other. The condition of the matrix can be improved 
and the matrix can be made symmetrical by the following transformation: 
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fi1 2 = 6 


(A12) 


which yields the final matrix 


- p‘a-lp) ^e" 1 



0 


(A13) 


The natural frequencies are then obtained when the determinant of 
equation (A13) vanishes. When the natural frequencies are substituted back 
into the matrix (eq. (A13)) and the value of unity assigned to one of 
the coefficients (for example, f^), the rest of the coefficients 

| 2 > and so forth, can be obtained. If an iteration procedure is used, 

the values of | are determined simultaneously with the natural fre- 
quencies. In either case, the mode shapes are then obtained as follows: 

n-L = -a -1 p ^ -1 i (Al^) 

and 

T12 = v/e" 1 ! (A15) 


The relative magnitudes of the elements of equation (A13) are illus- 
trated by the following matrix which results when ft = 30 radians per 
second and the parameters are those for the basic configuration: 


0.0001715 - -k 
of 

-0.00002270 

0.00008273 

-0.00001944 

0.0001195 

-0.00002808 


-0.00002270 
0.00004242 - i. 
-0.00000044-35 
0.0000001042 

-0.00003199 
0.00000751 8 


0.00008273 -0.00001944 

-0.0000004435 0.0000001042 
0.001030 - -i- 0.00001907 
0.00001907 O.OOOI563 - 

of 

0.00002841 -O.OOOOO6676 
-O.OOOOO6676 O.OOOOOI569 


0.0001195 -0.00002808 


f 

*1 

-0.00003199 0.000007518 


^2 

0.00002841 -O.OOOOO6676 


b 

-O.OOOOO6676 O.OOOOOI569 


*4 

0.001039 - ~ 0.00001684 


Ik 



5 

0.00001684 0.0001568 - i 


5 6 
L J 


(A16) 
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Procedure When a = 0 

When a - 0, the frequency parameter l/aP- occurs in only the last 
four rows and columns of equation (ll) , The transformation of the matrix 
for the solution in this case is accomplished in the same fashion as for 
U 4 0 except that the submatrices (a, (3, p t, e, and rjg of 

eq. (Al)) are of the fourth order. The relative magnitudes of the elements 
of the final matrix for Q = 0 and for the basic parameters are shown by 
the following equation: 




NACA TN 3849 


29 


REFERENCES 


1. Houbolt, John C. , and Anderson, Roger A.: Calculation of Uncoupled 

Modes and Frequencies in Bending or Torsion of Nonuniform Beams. 

NACA TN 1522, 1948. 

2. Myklestad, N. 0.: Vibration Analysis. McGraw-Hill Book Co., Inc., 

1944, pp. 184-214. 

3. Brooks, George W. : Analytical Determination of the Natural Coupled 

Frequencies of Tandem Helicopters. Jour. Am. Helicopter Soc., vol. 1, 
no. 3 , July 1956, pp. 39-52. 

4. Yntema, Robert T. : Simplified Procedures and Charts for the Rapid 

Estimation of Bending Frequencies of Rotating Beams. NACA TN 3459 ^ 
1955. (Supersedes NACA EM L54G02.) 



NACA TN 384-9 


TABLE I.- PARAMETERS USED IN CALCULATIONS 


(a) Blade parameters (blades on. front and rear rotor 
assumed to be identical) 

oip^i = = 25.3 radians/sec; Kp^ = K R ^]_ = 6.2; Rp = R R = 210 in.; 

P = 1= 3; Mp = M R = 0.194- lb-sea^/in.J 


Station 


Radial 

location, 

‘F or r R> 

in. 

Mass 

distribution, 
mp or m R , 

lb-sec^/in.^ 

First-bending 
mode shape, 

X 1 or z l> 
in. 

10.5 

3.03 

-0.09 

31.5 

1.52 

-.27 

52.5 

.67 

-.42 

73.5 

•59 

-.50 

94.5 

•59 

-•53 

115.5 

•53 

-.50 

136.5 

.45 

-.34 

157-5 

.40 

-.07 

178.5 

• 50 

.28 

199.5 

• 59 

•75 


(b) Fuselage and engine parameters 

-j. = 89.2 radians/sec; = 130 radians/sec; 
Me = 2.861 lb-sec2/in.; l e /l = O. 735 ] 


Fuselage 

location, 

s/l 


Length of Mass of First-bending 

fuselage element, fuselage element, mode shape. 
As, mf As, y x , 

lb-sec^/in. in. 


32.4-4- 

32.44 

32.44 

32.44 

32.44 

24.46 

40.42 

32.44 

32.44 


1.28 
1.98 . 
.42 
.86 
2.12 
.So 
• 23 
.81 
1.58 


O.76 

.36 

.03 

-.21 

-.36 

-.38 

-.27 

.26 

1.00 
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TABLE II.- MODAL COEFFICIENTS 



Modal coefficients for natural frequencies - 

®L 



• 


% * 

— ”1 


ft = 0 




-2.362 

86.36 

-0.02010 

0.02635 



- 0.6313 

- 157.3 

0.1598 

- 1.736 



2.652 

- 144.2 

- 1.224 

-1.301 



3-299 

115-7 

-1.886 

1.174 



1 

1 

1 

1 



-182.1 

2,802 

1.043 

1.032 



-226.5 

-2,250 

1.607 

- 0.9300 



- 0.1233 

-1.212 

-0.1182 

4.266 

ft = 20 radians/sec 

- 3-925 

123.1 

-0.2716 

32.21 

-0.01219 

0.03148 

-O.9383 

- 225.3 

- 0.04228 

-56.51 

0.1701 

- 1.748 

153.4 

-2,183 

- 0.9310 

-62.22 

- 1.310 

-1.330 

187.5 

1,782 

-1.308 

43.95 

-2.028 

1.204 

1 

1 

1 

1 

1 

1 

0.6192 

-26.11 

- 21.40 

1,060 

1.684 

1.148 

0.7558 

21.31 

-30.06 

43-95 

2.608 

- 1.040 

-O.1227 

- 1.146 

-0.1252 

-2.282 

-0.1171 

4.273 

ft = 30 radians /sec 

- 1.552 

62.72 

-O.O7493 

35-59 

0.009242 

0.04096 

- 0.4512 

-113.6 

0.09865 

-60.30 

0.2010 

- 1.775 

61.00 

- 1,123 

- 1.369 

-69.86 

- 1.425 

-1.366 

77.30 

881.8 

- 2.045 

45.57 

- 2.242 

1.266 

1 

1 

1 

1 

1 

l 

0.1760 

- 15-41 

- 4.654 

1,188 

3.734 

1.389 

0.2225 

12.11 

-6.968 

- 776.4 

5.883 

-1.274 

- 0.1239 

- 1.312 

- 0.1223 

-5.586 

- 0.1128 

4.293 
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Figure 2.- Variation of natural frequencies with rotor speed. 
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(d) Fourth— coupled— mode shape. 




I 


a 


L 


Figure 3-- Continued. 
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.(f) Sixths coupled- mode shape. 
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Figure 3«- Concluded. 
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Figure 5 .- Effect of Southwell coefficient for blade first elastic 
flapwise bending on natural frequencies of coupled system. 

0 = 30 radians per second. 
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Rotor speed, £1 , rodions/sec 

Figure 6.- Comparison of natural frequencies of coupled system from 
three-mode analysis with corresponding frequencies from eight-mode 
analysis. 
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o In-phase component, Xj 
□ Out-of-phase component, 


Natural coupled frequencies 


a>, radians/sec 


Figure 8.- Effect of forcing frequency on in-phase and out-of -phase 
components of amplitude of vibration of cockpit. Oscillating 
force of unit magnitude applied at rear rotor hub. 
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Figure 9 .- Effect of forcing frequency on amplitude of vibration of 
cockpit. Oscillating force of unit magnitude applied at front 
rotor hub. 


NACA TN 3849 


43 



Figure 10. - Effect of forcing frequency on amplitude of vibration of 
cockpit. Oscillating force of unit magnitude applied at rear rotor 
hub. 


<a, radians/sec 


Figure 11.- Effect of forcing frequency on amplitude of vibration of 
cockpit. Oscillating force of unit magnitude equally divided 
between front and rear rotor hubs and applied in same direction. 
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Figure 12 .- Effect of forcing frequency on amplitude of vibration of 
cockpit. Oscillating force of unit magnitude equally divided 
between front and rear rotor hubs and applied in opposite directions. 
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